Method of fault location in parallel lines with series compensation

ABSTRACT

The present invention relates to a method for locating a fault (F) in a section of parallel transmission lines in a network comprising the steps: measuring the currents and voltages of both lines at a measuring point arranged at one end (A) of the section, determining the fault distance (x) between the measuring point and the fault as a solution of an equation Ax 2 −Bx+C−R f =0 comprising the fault distance (x) as a variable and the fault resistance (R F ), the invention is characterized in that the parameters (A, B, C, D) comprise the phase components of the locally measured currents and voltages and are obtained from calculating from the measuring point to the fault location along the both parallel lines, and wherein the equation is resolved into its real and imaginary parts:  
     Real( A ) x   2 −Real( B ) x+ Real( C )− R   f =0  
     Imag( A ) x   2 −Imag( B ) x+ Imag( C )=0,  
     whereby the fault distance is derived from the imaginary part.

TECHNICAL FIELD

[0001] The present invention relates to a method for locating a fault (F) in a section of parallel transmission lines in a network comprising the steps, measuring the currents and voltages of both lines at a measuring point arranged at one end (A) of the section, determining the fault distance (x) between the measuring point and the fault as a solution of an equation comprising the fault distance as a variable and the fault resistance.

BACKGROUND OF THE INVENTION

[0002] Parallel series-compensated lines, i. e., lines provided with series capacitors and metal oxide varistors for improving the power transfer and enhancing power and voltage control of long transmission lines, are very important links between power generation and energy consumption regions. However, installation of series capacitors and varistors causes certain problems for the fault location.

[0003] Accurate fault location in parallel power transmission lines with series compensation requires compensating for the following effects:

[0004] 1. the remote infeed effect under resistive faults,

[0005] 2. the effect of the mutual coupling between the lines for the zero sequence,

[0006] 3. the effect of series compensation.

[0007] The countermeasures for the first two effects (1 and 2 as listed above) are applied for example in the single-ended fault locator proposed in U.S. Pat. No. 4,559,491. However, this fault locator is designated to locating faults in uncompensated lines, i. e., lines without series capacitor compensation. The cited method utilizes the local measurements (phase voltages, post-fault and pre-fault phase currents from the faulted line and a zero sequence current from the healthy line) as well as the impedance parameters for the lines and for the equivalent supplying systems at both the line ends. As the remote system impedance is not measurable with the single ended method, the fault locator according to U.S. Pat. No. 4,559,491 applies the representative value of this impedance for the positive sequence.

[0008] This is possible due to comparatively high robustness of this algorithm against mismatch of the actual and the representative values. However, in extreme cases of the high mismatch this is the additional source of error of the location algorithm. Moreover, the source impedances are also subjected to changes under evolving faults. In addition, if there is an extra link between the stations the impedance of this equivalent link ought to be provided for the location algorithm and obviously inaccuracy in this data affects the fault location too.

[0009] The countermeasure for the effect under point 3 above (the effect of series compensation) has been proposed in patent application PCT/SE98/02404, where the idea from U.S. Pat. No. 4,559,491, as well as in Article “A new fault locating algorithm for series compensated lines”, IEEE Trensactions on Power Delivery, Vol. 14, No. 3, July 1999, pp 789-795, is extended to the case of locating faults in a single line with series compensation. For this purpose the fundamental frequency equivalenting of parallel branches of a compensating series capacitor (SC) and a Metal Oxide Varistor (MOV) has been introduced. Currents flowing through MOVs in particular phases during unsymmetrical faults have different amplitudes. As a consequence of that the parameters (resistances and reactances) of the fundamental frequency equivalents are different in particular phases. So, MOVs in particular phases may have different fundamental frequency representations in the fault location algorithm.

BRIEF DESCRIPTION OF THE INVENTION

[0010] The aim of the present invention is to provide a fault location method which does not require the knowledge of the source impedances of the system behind both stations and takes into consideration the reactance effect, the series compensation effect and the mutual coupling between the lines.

[0011] This aim is obtained by the characterising part of claim 1.

[0012] The method according to the present invention displays advantages in relation to the above mentioned methods. It is suited for fault location in parallel lines with series compensation as well as after adequate setting for fault location in parallel uncompensated lines. It is based on phase coordinates approach which allows the incorporation of the fundamental frequency equivalents of the SCs&MOVs and to locate faults in the untransposed parallel lines. Further, it utilizes the local post-fault measurements, i.e., for the fault locator from the station A;—the phase voltages,—the phase currents from the faulted line and—the phase currents from the healthy line. It requires knowing the impedance parameters only for the lines (in terms of the self and mutual impedances of lines as well as for the mutual coupling between lines). Due to utilizing the healthy line path the impedances of the equivalent systems behind both the stations and the impedance of the equivalent link between the stations are not needed.

[0013] For determining the fault distance for series compensated lines, two subroutines are used: subroutine 1—estimating the distance to fault under assumption that it occurs behind the SCs&MOVs, and subroutine 2—estimating the distance to fault under assumption that it is applied in front of the SCs&MOVs. It is to be noted that usage of the algorithm for locating faults in parallel uncompensated lines relies in using only one of the subprocedures in which the fundamental frequency equivalents of the SCs&MOVs are set to zero for both the resistance and reactance parameters.

[0014] The sought fault distance is obtained with the extra procedure selecting the final result from the results of the both subroutines; the selecting procedure yields the indication which the subroutine is valid in a particular case on the base of information obtained:

[0015] a) the estimated fault resistances by both the subroutines,

[0016] b) amplitudes of the estimated currents of the healthy phases in the fault paths (for all the fault types except the three phase faults for which the selection procedure is performed only by (a)).

[0017] Thus, the proposed fault location method is an extension of the uncompensated lines fault locator described in U.S. Pat. No. 4,559,491 to the case of parallel lines with series compensation or can be treated as an extension of the single series compensated line fault locator according to PCT/SE98/02404 to the case of parallel series compensated lines. The main advantages are that it does not require the knowledge of the source impedances of the systems behind both the stations and the impedance of the equivalent link between the stations. Moreover, it does not use the pre-fault measurements.

[0018] The introduced phase coordinates approach makes possible to incorporate the fundamental frequency equivalents into the model as well as to consider a line as an untransposed line.

[0019] These and other aspects of, and advantages with the present invention, will become apparent form the detailed description of an embodiment and from the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In the detailed description reference will be made to the drawings, of which:

[0021]FIG. 1 is a schematic arrangement of series compensated parallel transmission lines,

[0022]FIG. 2a is the scheme of fundamental frequency equivalence of the series capacitors and metal-oxide varistors,

[0023]FIG. 2b shows a graph over the equivalent series resistance and reactance as a function of the amplitude of a current through the scheme of FIG. 2a,

[0024]FIG. 3 is a schematic arrangement for representing different types of faults,

[0025]FIG. 4 is a model of the system for faults located behind the series compensation as seen from the measuring point,

[0026]FIG. 5 is the model according to FIG. 4 where the fault is located in front of the series compensation,

[0027]FIG. 6 shows one example of the connection of a fault locator according to the invention to an existing line protection device, and

[0028]FIG. 7 shows an example of a device and system for carrying out the method.

DETAILED DESCRIPTION OF THE INVENTION

[0029] A model of the system for introducing the phase coordinates approach is presented in FIG. 1 it is to be noted that in this figure and along the whole this document the matrix and vector quantities are bold-typed.

[0030] All the quantities (voltages, currents) are represented further as vectors. For example, the side A phase voltages are represented by the vector: $\begin{matrix} {V_{A} = \begin{bmatrix} {\underset{\_}{V}}_{AR} \\ {\underset{\_}{V}}_{AS} \\ {\underset{\_}{V}}_{AT} \end{bmatrix}} & (1) \end{matrix}$

[0031] Impedance parameters are expressed in the matrix form. For example, the matrix of the impedance for a line takes the form: $\begin{matrix} {Z_{L} = \begin{bmatrix} {\underset{\_}{Z}}_{RR} & {\underset{\_}{Z}}_{RS} & {\underset{\_}{Z}}_{RT} \\ {\underset{\_}{Z}}_{SR} & {\underset{\_}{Z}}_{SS} & {\underset{\_}{Z}}_{ST} \\ {\underset{\_}{Z}}_{TR} & {\underset{\_}{Z}}_{TS} & {\underset{\_}{Z}}_{TT} \end{bmatrix}} & (2) \end{matrix}$

[0032] Self and mutual impedances from equation (2) representing completely symmetrical 3-phase line, are defined by the positive (Z_(L1)) and zero the (Z_(L0)) sequence impedances of a line as: $\begin{matrix} {{\underset{\_}{Z}}_{RR} = {{\underset{\_}{Z}}_{SS} = {{\underset{\_}{Z}}_{TT} = \frac{{\underset{\_}{Z}}_{L0} + {\underset{\_}{Z}}_{L1}}{3}}}} & (3) \\ {{\underset{\_}{Z}}_{RS} = {{\underset{\_}{Z}}_{RT} = {{\underset{\_}{Z}}_{SR} = {{\underset{\_}{Z}}_{ST} = {{\underset{\_}{Z}}_{TR} = {{\underset{\_}{Z}}_{TS} = \frac{{\underset{\_}{Z}}_{L0} - {\underset{\_}{Z}}_{L1}}{3}}}}}}} & (4) \end{matrix}$

[0033] The SC&MOV circuit in each phase of a line is linearized for the steady-state conditions and represented in the form of the fundamental frequency equivalent series R_(v)−X_(v) connection (FIG. 2a).

[0034] The matrix of equivalent parameters for the SCs&MOVs from the line A (the current I_(AA)) is given as: $\begin{matrix} {{Z_{v}\left( {{\underset{\_}{I}}_{AA}} \right)} = \begin{bmatrix} {{\underset{\_}{Z}}_{v}\left( {{\underset{\_}{I}}_{AA\_ R}} \right)} & 0 & 0 \\ 0 & {{\underset{\_}{Z}}_{v}\left( {{\underset{\_}{I}}_{AA\_ S}} \right)} & 0 \\ 0 & 0 & {{\underset{\_}{Z}}_{v}\left( {{\underset{\_}{I}}_{AA\_ T}} \right)} \end{bmatrix}} & (5) \end{matrix}$

[0035] Diagonal elements of the matrix Z_(v) depend on magnitudes of phase currents. Their real (resistance: R_(v)) and imaginary (reactance: X_(v)) components are obtained from the relations graphically shown in FIG. 2b.

[0036] Faults in the considered system are described using the 3-phase general fault model (FIG. 3), which is representing a general fault situation. For this scheme one gets:

G _(f) V _(f) =I _(f)  (6)

[0037] where: $G_{f} = \begin{bmatrix} G_{RR} & {- G_{RS}} & {- G_{RT}} \\ {- G_{RS}} & G_{SS} & {- G_{ST}} \\ {- G_{RT}} & {- G_{ST}} & G_{TT} \end{bmatrix}$ $G_{RR} = {\frac{1}{R_{R}} + \frac{1}{R_{RS}} + \frac{1}{R_{RT}}}$ $\begin{matrix} {G_{RS} = \frac{1}{R_{RS}}} & {G_{RT} = \frac{1}{R_{RT}}} \\ {V_{f} = \begin{bmatrix} {\underset{\_}{V}}_{fR} \\ {\underset{\_}{V}}_{fS} \\ {\underset{\_}{V}}_{fT} \end{bmatrix}} & {I_{f} = \begin{bmatrix} {\underset{\_}{I}}_{fR} \\ {\underset{\_}{I}}_{fS} \\ {\underset{\_}{I}}_{fT} \end{bmatrix}} \end{matrix}$

[0038] Assuming all the resistances in the fault circuit to be the same, say R_(f), matrix G_(f) takes the following form: $\begin{matrix} {{G_{f} = {\frac{1}{R_{f}}K_{f}}}{K_{f} = \begin{bmatrix} k_{RR} & k_{RS} & k_{RT} \\ k_{RS} & k_{SS} & k_{ST} \\ k_{RT} & k_{ST} & k_{TT} \end{bmatrix}}} & (7) \end{matrix}$

[0039] The elements in matrix K_(f) are determined in dependence on the type of fault as follows:

[0040] non-diagonal elements are given the value 0 if the phase in question is not concerned by the relevant fault and the value −1 if the phase is concerned by the relevant fault

[0041] the diagonal elements are given the value 1 if the phase in question has a fault to ground at the fault in question and to this added the sum of the absolute values of the non-diagonal elements in the relevant line.

[0042] A few examples of a filled-in matrix K_(f) for some typical types of fault are shown below. $K_{f} = {\begin{bmatrix} 1 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix}\quad {for}\quad {fault}\quad {to}\quad {ground}}$ $K_{f} = {\begin{bmatrix} 1 & {- 1} & 0 \\ {- 1} & 1 & 0 \\ 0 & 0 & 0 \end{bmatrix}\quad {for}\quad {fault}\quad R\text{-}S}$ $K_{f} = {\left\lbrack \quad \begin{matrix} 3 & {- 1} & {- 1} \\ {- 1} & 3 & {- 1} \\ {- 1} & {- 1} & 3 \end{matrix}\quad \right\rbrack \quad {for}\quad {fault}\quad \text{R-S-T-G}}$ $K_{f} = {\begin{bmatrix} 2 & {- 1} & 0 \\ {- 1} & 2 & 0 \\ 0 & 0 & 0 \end{bmatrix}\quad {for}\quad {fault}\quad \text{R-S-G}}$ $K_{f} = {\begin{bmatrix} 2 & {- 1} & {- 1} \\ {- 1} & 2 & {- 1} \\ {- 1} & {- 1} & 2 \end{bmatrix}\quad {for}\quad {fault}\quad \text{R-S-T}}$

[0043] Fault Location Algorithm

[0044] Below, a fault location algorithm is derived assuming that fault location process is performed with use of the post-fault quantities from the station A and for the faulted line A (FIGS. 4,5). There are two distinctive fault locations considered by the fault location algorithm with the subroutines considering:

[0045] faults behind the SCs&MOVs (Subroutine 1)

[0046] faults in front of the SCs&MOVs (Subroutine 2)

[0047] Fault Location Algorithm—the Subroutine for the Case of a Fault Behind SCs&MOVs (SUBROUTINE 1)

[0048]FIG. 4 shows the series compensated network for the case of faults occurring behind the SCs&MOVs, but not overreaching the total line length. In this case a fault loop as seen from the station A contains the SCs&MOVs and the infeed is via the remote segment of the faulted line. In addition, mutual coupling between the healthy and the faulted lines has to be taken into account.

[0049] Using matrix notation, the model of the series-compensated network from FIG. 4 is described:

V _(f) =V _(A) −xZ _(LA) I _(AA) −Z _(v)(|I _(AA)|)I _(AA) −xZ _(m) I _(AB)  (8a)

[0050] or:

V _(f) =V _(B)−(1x)Z _(LA) I _(BA+()1−x)Z _(m) I _(AB)  (8b)

[0051] by utilizing:

I _(BB) =−I _(AB)

[0052] where:

[0053] x—unknown and the sought distance to a fault, counted from the station A to the fault place

(p≦X≦1),

[0054] Z_(LA)—impedance matrix for the faulted line,

[0055] Z_(v)(|I_(AA)|)—three-phase fundamental frequency equivalent of the SCs&MOVs from the faulted line,

[0056] Z_(m)—impedance matrix for the mutual coupling between the lines.

[0057] It is to be noted that for the healthy line, equation (8b), it is assumed that phase currents at both the line ends are equal to each other (this is a consequence of neglecting shunt capacitances).

[0058] On the other hand, the voltage drop between the stations A and B can be calculated by considering the path of the healthy line:

V _(A) −V _(B)=(Z _(LB) +Z _(v)(|I _(AB)|))I _(AB) +xZ _(m) I _(AA)−(1−x)Z _(m) I _(BA)  (8c)

[0059] where:

[0060] Z_(LB)−impedance matrix for the healthy line (basically equal to the impedance matrix of the faulty line Z_(LA)),

[0061] Z_(v)(|I_(AB)|)—three-phase fundamental frequency equivalent of the SCs&MOVs from the healthy line.

[0062] Subtracting equation (8a) and equation (8b) gives

0=( V _(A) −V _(B))−xZ _(1A) I _(AA) −Z _(v)(|I _(AA)|)I _(AA) +Z _(1A) I _(BA) −Z _(m) I _(AB)  (9)

[0063] Substituting equation (8c) into equation (9) and rearranging gives three-phase currents flowing in the faulted line in the remote station B: $\begin{matrix} {I_{BA} = {A_{0} + \frac{{xB}_{0}I_{AA}}{B_{0}\left( {1 - x} \right)}}} & (10) \end{matrix}$

[0064] where:

[0065] A₀=(Z_(LB)+Z_(v)(|I_(AB)|)−Z_(m))I_(AB)−Z_(v)(|I_(AA)|)I_(AA)

[0066] B₀=Z_(m)−Z_(LA)

[0067] Voltages at a fault point (FIG. 4) can be expressed as in equation (8a):

V _(f) =V _(A) −xZ _(LA) I _(AA) −Z _(V)(|I _(AA)|)I _(AA) −xZ _(m) I _(AB)  (11)

[0068] If shunt capacitances of the line are neglected, which has been assumed here, the currents in the fault paths (FIG. 4) are expressed by a sum of the phase currents from both the stations:

I _(f) =I _(AA) +I _(BA)  (12)

[0069] Substituting equation (11) and equation (12) into equation (6) gives:

G _(f)(V _(A)−(xZ _(LA) +Z _(v)(|I _(AA)|))I _(AA) −xZ _(m) I _(AB))=I _(AA) +I _(BA)  (13)

[0070] Substituting remote station phase currents I_(BA) equation (10) into equation (13) and taking also equation (7), after performing multiplications and further rearranging, the following matrix formula (second order equation with respect to the sought distance to a fault, x) is obtained:

A _(c) x ² −B _(c) x+C _(c) −D _(c)=0  (14)

[0071] where:

[0072] A_(c), B_(c), C_(c), D_(c)−3*1 vectors:

[0073] A_(c)=(Z_(m)−Z_(LA))K_(f)(Z_(LA)I_(AA)+Z_(m)I_(AB))

[0074] C_(c)=(Z_(m)−Z_(LA))K_(f)(V_(A)−Z_(v)(|I_(AA)|)I_(AA))

[0075] B_(c)=A_(c)+C_(c)

[0076] D=(Z_(m)−Z_(LA)−Z_(v)(|I_(AA)|))I_(AA)−(Z_(m)−Z_(LB)−Z_(v)(|I_(AB)|))I_(AB)

[0077] D_(c)=DR_(f)

[0078] where:

[0079] R_(f)—equivalent fault resistance equation (7).

[0080] Equation (14) represents a matrix formula for the sought distance to a fault, x, and the unknown equivalent fault resistance R_(f). Multiplying both the sides of equation (14) by the vector: $\begin{matrix} {P = \frac{D}{D^{T}D}} & (15) \end{matrix}$

[0081] where:

[0082] D^(T)—matrix transposed with respect to the matrix D, one obtains the following resultant complex scalar equation:

Ax ² −Bx+C−R _(f)=0  (16)

[0083] where:

[0084] A=PA_(c)

[0085] B=PB_(c)

[0086] C=PC_(c)

[0087] P=D^(T)/(D^(T)D) (vector)

[0088] D^(T)—matrix transposed with respect to D,

[0089] The scalar quadratic equation (16) can be resolved into its real and imaginary parts:

Real(A)x ²−Real(B)x+Real(C)−R _(f)=0  (17)

Imag(A)x ²−Imag(B)x+Imag(C)−R _(f)=0  (18)

[0090] The imaginary part, equation (18), does not contain R_(f). Solving it, one obtains the sought distance to a fault, x₁ [pu], according to the SUBROUTINE 1, as:

x ₁ =x _(a), if Imag(A)>0,

x ₁ =x _(b), if Imag(A)<0,  (19)

[0091] where: $x_{a} = \frac{{{Imag}(B)} - \sqrt{D}}{2{{Imag}(A)}}$ $x_{b} = \frac{{{Imag}(B)} + \sqrt{D}}{2{{Imag}(A)}}$ Δ = [Imag(B)]² − 4Imag(A)Imag(C)

[0092] Knowing the distance to a fault, equation (19), (x₁=x_(a) or x₁=x_(b)) one can determine the value of the equivalent fault resistance (R_(f)) from equation (17) as:

R _(f)=Real(A)x ₁ ²−Real(B)x ₁+Real(C)  (20)

[0093] Fault Location Algorithm—the Subroutine for the Case of a Fault in Front of SCs&MOVs (SUBROUTINE 2)

[0094]FIG. 5 presents the series compensated network for the case of faults occurring in front of the SCs&MOVs. In this case a fault loop, as seen from the station A, does not contain the SCs&MOVs. The infeed is via the remote segment of the faulted line together with the SCs&MOVs. It is worth to notice that the SCs&MOVs even are not present in the fault loops considered by the fault location algorithm they influence the fault location. This is so, because the SC&MOVs influence the remote infeed. Therefore, the SCs&MOVs have to be taken into account when deriving the fault location subroutine of this case too.

[0095] As in the previous subroutine, using matrix notation, one can write the following equation for the faulted network of FIG. 5:

V _(f) =V _(A) −xZ _(LA) I _(AA) −xZ _(m) I _(AB)  (22a)

[0096] or:

V _(f) =V _(B)−(1−x)Z _(LA) I _(BA)+(1−x)Z _(m) I _(AB) −Z _(v)(|I _(BA)|)I _(BA)  (22b)

[0097] where:

[0098] x—unknown and the sought distance to a fault, counted from the station A to the fault place (0≦x≦p),

[0099] Z_(LA)—impedance matrix for the faulted line,

[0100] Z_(v)(|I_(BA)|)—three-phase fundamental frequency equivalent of the SCs&MOVs from the faulted line,

[0101] Z_(m)—impedance matrix for the mutual coupling between the lines.

[0102] It is to be noted that for the healthy line it is assumed that phase currents at both the line ends are equal to each other (this is a consequence of neglecting shunt capacitances).

[0103] On the other hand, the voltage drop between the stations A and B can be calculated by considering the path of the healthy line:

V _(A) −V _(B)=(Z _(LB) +Z _(v)(|I _(AB)|))I _(AB) +xZ _(m) I _(AA)−(1−x)Z _(m) I _(BA)  (23)

[0104] where:

[0105] Z_(LB)—impedance matrix for the healthy line,

[0106] Z_(v)(|I_(AB)|)—three-phase fundamental frequency equivalent of the SCs&MOVs from the healthy line.

[0107] Three-phase currents flowing in the faulted line at the remote station B can be determined from equations (22) and (23) as: $\begin{matrix} {I_{BA} = \frac{C_{0} + {{xB}_{0}I_{AA}}}{\left( {B_{0} - {Z_{v}\left( {I_{BA}} \right)}} \right) - {xB}_{0}}} & (24) \end{matrix}$

[0108] where:

C ₀=(Z _(LB) +Z _(v)(|I _(AB)|)−Z _(m))I _(AB)

B ₀ =Z _(m) −Z _(LA)

[0109] Voltages at a fault point (FIG. 5) can be expressed as in equation (22a):

V _(f) =V _(A) −xZ _(LA) I _(AA) −xZ _(m) I _(AB)  (25)

[0110] If shunt capacitances of the line are neglected (which has been assumed here) currents in the fault paths (FIG. 5) are expressed by a sum of the phase currents from both the stations:

I _(f) =I _(AA) +I _(BA)  (26)

[0111] Substituting equation (25) and equation (26) into equation (6) gives:

G _(f)(V _(A)−(xZ _(LA) +Z _(v)(|I _(AA)|))I _(AA) −xZ _(m) I _(AB))=I _(AA) +I _(BA)  (27)

[0112] Substituting the remote station phase currents I_(BA) in equation (24) into equation (27) and taking also equation (7) after performing multiplications and further rearranging, the following matrix formula (second order equation with respect to the sought distance to fault, x) is obtained:

A _(c) x ² −B _(c) x+C _(c) −D _(c)=0  (28)

[0113] where:

[0114] A_(c), B_(c), C_(c), D_(c)—3*1 vectors:

[0115] A_(c)=(Z_(m)−Z_(LA))K_(f)(Z_(LA)I_(AA)+Z_(m)I_(AB))

[0116] B_(c)=(Z_(m)−Z_(LA))K_(f)(V_(A)+Z_(LA)I_(AA)+Z_(m)I_(AB))−Z_(v)(|I_(BA)|)K_(f)(Z_(LA)I_(AA)+Z_(m)I_(AB))

[0117] C_(c)=(Z_(m)−Z_(LA)−Z_(v)(|I_(BA |)K) _(f)V_(A)

[0118] D=(Z_(m)−Z_(LA)−Z_(v)(|I_(BA)|))I_(AA)−(Z_(m)−Z_(LB −Z) _(v)(|I_(AB)|))I_(AB)

[0119] D_(c)=DR_(f)

[0120] where:

[0121] R_(f)—equivalent fault resistance, equation (7).

[0122] Equation (28) represents a matrix formula for the sought distance to a fault, x, and the unknown equivalent fault resistance R_(f). Multiplying both the sides of equation (28) by the vector: $\begin{matrix} {P = \frac{D}{D^{T}D}} & (29) \end{matrix}$

[0123] where:

[0124] D^(T)—matrix transposed with respect to the matrix D, one obtains the following resultant complex scalar equation:

Ax ² −Bx+C−R _(f)=0  (30)

[0125] where:

[0126] A=PA_(c)

[0127] B=PB_(c)

[0128] C=PC_(c) $P = \frac{D}{D^{T}D}$

[0129] The scalar quadratic equation (30) can be resolved into its real and imaginary parts:

Real(A)x ²−Real(B)x+Real(C)−R_(f)=0  (31)

Imag(A)x ²−Imag(B)x+Imag(C)=0  (32)

[0130] The imaginary part, equation (32), does not contain (R_(f)): Solving it, one obtains the sought distance to a fault, x₂ [pu], according to the SUBROUTINE 2, as:

x ₂ =x _(a), if Imag(A)>0,

x ₂ =x _(b), if Imag(A)<0,  (33)

[0131] where: $x_{a} = \frac{{{Imag}(B)} - \sqrt{D}}{2{{Imag}(A)}}$ $x_{b} = \frac{{{Imag}(B)} + \sqrt{D}}{2{{Imag}(A)}}$ Δ = [Imag(B)]² − 4Imag(A)Imag(C)

[0132] Knowing the distance to a fault, equation (33), (x₂=x_(a) or x₂=x_(b)) one can determine the equivalent fault resistance (R_(f)) from equation (31) as:

R_(f)=Real(A)x ₁ ²−Real(B)x ₁+Real(C)  (34)

[0133] Iterative Calculation of the Remote Phase Currents for the Subroutine Assuming a Fault in Front of the SCs&MOVs (Subroutine 2)

[0134] Solution of the resultant quadratic scalar equation (33) (for the SUBROUTINE 2

[0135] a fault assumed as occurring in front of the SCs&MOVs) requires knowing fundamental frequency three-phase equivalents of the SCs&MOVs. In case of the faulted line the amplitudes of the remote currents have to be known. This can be be achieved by means:

[0136] amplitudes of the remote phase currents are iteratively estimated (see next section),

[0137] amplitudes of the remote phase currents are sent via communication channel (note: there is no need for synchronization of the data collection at both the ends).

[0138] Iterative Algorithm for Estimation of the Amplitudes of the Remote Currents

[0139] Determination of the amplitudes of the remote currents can be performed with the iterative algorithm for solution of the following set of nonlinear equations (note: numbering of equations is kept as in the section 3.2 where the SUBROUTINE 2 is derived):

[0140] Equation for the remote currents $\begin{matrix} {I_{BA} = \frac{C_{0} + {{xB}_{0}I_{AA}}}{\left( {B_{0} - {Z_{v}\left( {I_{BA}} \right)}} \right) - {xB}_{0}}} & (24) \end{matrix}$

[0141] where:

[0142] C₀=(Z_(LB)+Z_(v)(|I_(AB)|)−Z_(m))I_(AB)

[0143] B₀=Z_(m)−Z_(LA)

[0144] Equation for the distance to a fault:

Imag(A)x ²−Imag(B)x+Imag(C)=0  (32)

[0145] The sought distance to a fault, x₂ [pu], according to the SUBROUTINE 2, as:

x ₂ =x _(a), if Imag(A)>0,

x ₂ =x _(b), if Imag(A)<0,  (33)

[0146] where: $x_{a} = \frac{{{Imag}(B)} - \sqrt{\Delta}}{2{{Imag}(A)}}$ $x_{b} = \frac{{{Imag}(B)} + \sqrt{\Delta}}{2{{Imag}(A)}}$ Δ = [Imag(B)]² − 4Imag(A)Imag(C)

[0147] A=PA_(c)

[0148] B=PB_(c)

[0149] C=PC_(c)

[0150] P=D^(T)/(D^(T)D) (vector)

[0151] D^(T)—matrix transposed with respect to D,

[0152] where:

[0153] A_(c), B_(c), C_(c), D_(c) —3*1 vectors:

[0154] A_(c)=(Z_(m)−Z_(LA))K_(f)(Z_(LA)I_(AA)+Z_(m)I_(AB))

[0155] B_(c)=(Z_(m)−Z_(LA))K_(f)(V_(A)+Z_(LA)I_(AA)+Z_(m)I_(AB))−Z_(v)(|I_(BA)|)K_(f)(Z_(LA)I_(AA)+Z_(m)I_(AB))

[0156] C_(c)=(Z_(m)−Z_(LA)−Z_(v)(|I_(BA)|)K_(f)V_(A)

[0157] D=(Z_(m)−Z_(LA)−Z_(v)(|I_(BA)|))I_(AA)−(Z_(m)−Z_(LB)−Z_(v)(|I_(AB)|))I_(AB)

[0158] Equation determining the fundamental frequency equivalent of the SCs&MOVs: $\begin{matrix} {{{\underset{\_}{Z}}_{v}\left( {{\underset{\_}{I}}_{BA}} \right)} = \begin{bmatrix} {{\underset{\_}{Z}}_{v}\left( {{\underset{\_}{I}}_{BA\_ R}} \right)} & 0 & 0 \\ 0 & {{\underset{\_}{Z}}_{v}\left( {{\underset{\_}{I}}_{BA\_ S}} \right)} & 0 \\ 0 & 0 & {{\underset{\_}{Z}}_{v}\left( {{\underset{\_}{I}}_{BA\_ T}} \right)} \end{bmatrix}} & (35) \end{matrix}$

[0159] Determining the distance to the fault is based in an iterative process. before the iteration begins, the impedance of the series compensation of the faulty line is set and the initial value of I_(BA) is calculated.

[0160] During the iteration process the coefficients A, B and C of the equation for fault distance calculation Ax²−Bx+C−R_(f)=0 are calculated and subsequently the distance to fault. The first calculated value of the distance to the fault is compared with the assumed value of the distance to the fault. If the difference between the values is greater than a pre-set maximally permissible difference value, the assumed value of I_(BA) is replaced. Thereafter a new calculation procedure is carried out according to above which provides new values. This iteration continues until the difference value of two consecutive calculation procedures is smaller that the set maximally permissible difference value, whereupon the last calculated values of the distance to fault is regarded as the real value.

[0161] Selection Between the Subroutines

[0162] A selection between the subprocedures is performed with the properly weighted function of the two quantities:

[0163] estimated fault resistances by the SUBPROCEDURE 1 and 2: R_(F—SUB1), R_(F—SUB2)

[0164] estimated amplitudes of the healthy phases fault currents by the SUBPROCEDURE 1 and 2: |I_(f—healthy—ph)|_(SUB1), |I_(f—healthy—ph)|_(SUB2); for the considered cases: the average from the phases S, T amplitudes for the R-G faults and the phase T amplitude for the R-S-G faults are taken into consideration.

[0165] Lower value (if positive) of the estimated fault resistance and lower amplitude of the estimated fault path current support the selection of a particular subprocedure.

[0166] A fault locator according to the present invention may be connected directly to the line section as shown in FIG. 1 or via a line protection device LS as shown in FIG. 6. The device and the fault locator are supplied with the measured values of the currents of the both lines I_(AA), I_(AB) and voltages V_(A) via voltage and current transformers (not shown). In the case of a fault appearing on line LA the line protecting device LS delivers a signal S to the fault locator. On the one hand this signal triggers the start of the calculation, by the fault locator, of the distance to the fault, and on the other hand the signal S contains information about the type of fault that has arisen. Where necessary the line protecting device may also provide a signal that trips a circuit breaker (not shown). If the fault locator is connected directly, it may be provided with its own triggering means.

[0167]FIG. 7 shows an embodiment of a device for determining the distance from a station, at one end of a transmission line, until the occurrence of a fault on the transmission line according to the described method, comprising certain measuring devices, measurement value converters, members for treatment of the calculating algorithms of the method, indicating means for the calculated distance to fault and a printer for printout of the calculated fault.

[0168] In the embodiment shown, measuring devices 1 to 3 for continuous measurement of all the phase currents for both the faulted line LA and healthy line LB and phase voltages are arranged in station A. In the measurement converters 4 to 6, a number of these consecutively measured values, which in case of a fault are passed to a calculating unit 7, are filtered and stored. The calculating unit is provided with the calculating algorithms described, programmed for the processes needed for calculating the distance to fault and the fault resistance. The calculating unit is also provided with known values such as the impedance of the lines. In connection to the occurrence of a fault information regarding the type of fault may be supplied to the calculating unit for choosing the right coefficients. When the calculating unit has determined the distance to fault, it is displayed on the device and/or sent to remotely located display means. A printout of the result may also be provided. In addition to signalling the fault distance, the device can produce reports in which are recorded measured values of the currents of both lines, voltages, type of fault and other associated with a given fault at a distance.

[0169] The calculating unit may comprise filters for filtering the signals, A/D-converters for converting and sampling the signals and one or more micro processors. The micro processor (or processors) comprises a central processing unit CPU performing the following functions: collection of measured values, processing of measured values, calculation of distance to fault and output of result from calculation. The micro processor (or processors) further comprises a data memory and a program memory.

[0170] A computer program for carrying out the method according to the present invention is stored in the program memory. It is to be understood that the computer program may also be run on general purpose computer instead of a specially adapted computer.

[0171] The software includes computer program code elements or software code portions that make the computer perform the said method using the equations, algorithms, data and calculations previously described. A part of the program may be stored in a processor as above, but also in a RAM, ROM, PROM or EPROM chip or similar. The program in part or in whole may also be stored on, or in, other suitable computer readable medium such as a magnetic disk, CD-ROM or DVD disk, hard disk, magneto-optical memory storage means, in volatile memory, in flash memory, as firmware, or stored on a data server.

[0172] It is to be understood that the embodiments described above and shown on the drawings are to be regarded as non-limiting examples of the present invention and that it is defined by the appended patent claims. 

1. Method for locating a fault (F) in a section of parallel transmission lines in a network comprising the steps: measuring the currents and voltages of both lines at a measuring point arranged at one end (A) of the section, determining the fault distance (x) between the measuring point and the fault as a solution of an equation comprising the fault distance (x) as a variable and the fault resistance (R_(F)), characterised in that the equation is Ax ² −Bx+C−R _(f)=0 wherein the parameters (A, B, C, D) comprise the phase components of the locally measured currents and voltages and are obtained from calculating from the measuring point to the fault location along the both parallel lines, and wherein the equation is resolved into its real and imaginary parts: Real(A)x²−Real(B)x+Real(C)−R_(f)0 Imag(A)x ²−Imag(B)x+Imag(C)=0, whereby the fault distance is derived from the imaginary part, as x=x _(a), if Imag(A)>0, x=x _(b), if Imag(A)<0, where: $x_{a} = \frac{{{Imag}(B)} - \sqrt{D}}{2{{Imag}(A)}}$ $x_{b} = \frac{{{Imag}(B)} + \sqrt{D}}{2{{Imag}(A)}}$ Δ = [Imag(B)]² − 4Imag(A)Imag(C).


2. Method according to claim 1, characterised in that the parameters comprise the particular type of fault.
 3. Method according to claim 2, characterised in determining a matrix K_(f) for the particular type of fault by using a 3-phase general fault model.
 4. Method according to claim 3, characterised in the further steps: using a matrix notation of the section taking into account of the mutual impedances between the lines, including the fault type matrix, whereby obtaining a matrix formula A _(c) x ² −B _(c) x+C _(c) −D _(c)=0 where A_(c), B_(c), C_(c), D_(c) are 3*1 vectors, and multiplying both sides of the matrix formula with vector $P = \frac{D}{D^{T}D}$

where D^(T) is matrix transposed with respect to matrix D and A_(c)=(Z_(m)Z_(LA))K_(f)(Z_(LA)I_(AA)+Z_(m)I_(AB)) C_(c)=(Z_(m)−Z_(LA))K_(f)(V_(A)−Z_(v)(|I_(AA)|)I_(AA)) B_(c)=A_(c)+C_(c) D=(Z_(m)−Z_(LA)−Z_(v)(|I_(AA)|))I_(AA)−(Z_(m)−Z_(LB)−Z_(v)(|I_(AB)|))I_(AB) D_(c)=DR_(f)
 5. Method according to any of the preceding claims, characterised in that the matrix K_(f) is expressed as $K_{f} = \begin{bmatrix} k_{RR} & k_{RS} & k_{RT} \\ k_{RS} & k_{SS} & k_{ST} \\ k_{RT} & k_{ST} & k_{TT} \end{bmatrix}$


6. Method according to any of the preceding claims, characterised in that the parallel transmission lines are series compensated and that the compensation in the calculating paths is represented as equivalent resistance and reactance.
 7. Method according to claim 6, characterised in that the parallel connection of a series capacitor and a varistor constitutes a non-linear impedance (Z_(V)) which is represented by the equivalent resistance and reactance (R_(V) and X_(V)) which are determined as a function of a traversing current of each phase, whereupon the actual value of the resistance and reactance may be determined with the actual currents, which, after the occurrence of the fault, flow through the impedance, whereby the non-linear impedance may be set in a matrix form.
 8. Method according to claim 6 or 7, characterised in solving the equation for two cases; (1) when the fault is assumed behind the series compensation of the faulted line and (2) when the fault is assumed in front of the series compensation of the faulted line, as seen from the measuring point.
 9. Method according to claim 6 or 7, characterised in determining a matrix of equivalent parameters for the series capacitors and movistors of the faulted line in order to obtain the fundamental frequency equivalent.
 10. Method according to claim 9, characterised in calculating the fault resistance R_(f) from the real part of the equation and including the matrix of the fundamental frequency.
 11. Method according to claim 10, characterised in that the matrix of the fundamental frequency is expressed as ${Z_{v}\left( {I_{AA}} \right)} = \begin{bmatrix} {{\underset{\_}{Z}}_{v}\left( {I_{{AA}\quad \_ \quad R}} \right)} & 0 & 0 \\ 0 & {{\underset{\_}{Z}}_{v}\left( {I_{{AA}\quad \_ \quad S}} \right)} & 0 \\ 0 & 0 & {{\underset{\_}{Z}}_{v}\left( {I_{{AA}\quad \_ \quad T}} \right)} \end{bmatrix}$

for case (1), and ${Z_{v}\left( {I_{BA}} \right)} = \begin{bmatrix} {{\underset{\_}{Z}}_{v}\left( {I_{{BA}\quad \_ \quad R}} \right)} & 0 & 0 \\ 0 & {{\underset{\_}{Z}}_{v}\left( {I_{{BA}\quad \_ \quad S}} \right)} & 0 \\ 0 & 0 & {{\underset{\_}{Z}}_{v}\left( {I_{{BA}\quad \_ \quad T}} \right)} \end{bmatrix}$

for case (2).
 12. Method according to claim 10, characterised in selecting the calculated fault distance from the two cases, based on estimated fault resistance of the two cases, and estimated amplitudes of the healthy phases fault current, whereby lower estimated values support the calculated fault distance of a case.
 13. Method according to claim 12, characterised in that the amplitudes of the phase currents for case (2) are obtained by iterative calculation.
 14. Device for locating a fault (F) in a section of parallel transmission lines comprising calculating members arranged to calculate, on the basis of current and voltage values measured adjacent to one end of said section and the known impedance of the lines, the distance between the measuring point and the fault, characterised in that the calculating members are arranged to determine, on the basis of information about the type of fault in question and the complex quantities of the measured values and using a network model, the fault distance as the solution of an equation Ax ² −Bx+C−R _(f)=0 wherein the parameters (A, B, C, D) comprise the phase components of the locally measured currents and voltages and are obtained from calculating from the measuring point to the fault location along the both parallel lines, and wherein the equation is resolved into its real and imaginary parts: Real(A)x ²−Real(B)x+Real(C)−R _(f)=0 Imag(A)x ²−Imag(B)x+Imag(C)=0, whereby the fault distance is derived from the imaginary part, as x=x _(a), if Imag(A)>0, x=x _(b), if Imag(A)<0, where: $x_{a} = \frac{{{Imag}(B)} - \sqrt{D}}{2{{Imag}(A)}}$ $x_{b} = \frac{{{Imag}(B)} + \sqrt{D}}{2{{Imag}(A)}}$ Δ   = [Imag(B)]² − 4Imag(A)Imag(C).


15. Use of a device according to claim 14 to determine the distance to fault in a parallel transmission line.
 16. Use of a device according to claim 14 to record and signal currents and voltages associated with a fault at a distance (d) from a measuring station (A).
 17. Computer program product comprising computer code means and/or software code portions for making a computer or processor perform the steps of: receiving values of currents and voltages of both the lines at a measuring point arranged at one end (A) of the section, calculating a fault distance (d) from the measuring point to the fault location along the both parallel lines as a solution of an equation comprising the fault distance (d) as a variable and the fault resistance (Rf) and parameters that depend on the phase components of the locally measured currents and voltages, wherein the equation is resolved for the real parts and the imaginary parts respectively, and reporting the fault distance (d).
 18. Computer program product according to claim 17 contained on, or in, a computer readable medium. 